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STATEMENT OF SCOPE AND PURPOSE 



This paper presents new computational methods that facilitate digital com- 
puter analysis of some important military operations research problems. Even though 
combat between two military forces is a complex random process, deterministic 
Lanchester-type differential equations [1] are commonly used in defense planning 
studies. A so-called attrition-rate coefficient in such a combat model represents 
the fire effectiveness of a weapon-system type against a particular target type, 
i.e. its effective firepower. Time-dependent attrition-rate coefficients are used 
to model temporal variations in firepower on the battlefield. For such a variable- 
coefficient Lanchester-type combat model, we present a simple numerical procedure 
that allows one (without having to explicitly solve the equations) not only to predict 
battle outcome but also to parametrically tradeoff quality versus quantity of two 
opposing weapon systems. 
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ABSTRACT 



This paper presents a simple numerical procedure for determining the 
parity-condition parameter for Lanches ter- type combat between two homogeneous 
forces. The combat studied is modelled by Lanches ter- type equations of modern 
warfare with time-dependent attrition-rate coefficients. Previous research has 
shown that the prediction of battle outcome (in particular, force annihilation) 
without having to spend the time and effort of computing force-level trajectories 
depends on a single parameter, the so-called parity-condition parameter, which 
only depends on the attrition-rate coefficients. Unfortunately, previous research 
did not show how to generally determine this parameter. We present general 
theoretical considerations for its numerical noniterative determination. This 
general theory is applied to an important class of attrition-rate coefficients 
(offset power attrition-rate coefficients) . Our results allow one to study 
such variable-coefficient combat models almost as easily and thoroughly as 
Lanches ter *s classic constant-coefficient model. 
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1. Introduction 



As a consequence of pioneering work by F. W. Lanchester [16] done about the time 

• 4 - 

of World War I, military operations analysts have used simplified deterministic dif- 
ferential-equation models to develop insights into the dynamics of combat from about 
the end of World War II (see, for example, [1], [7-8], [13], [24-28]). The advent of 
the modern high-speed digital computer has made feasible the development and use of 
quite complicated versions of such Lanchester-type^ models as practical defense planning 
tools [6]. Thus, today militarily realistic computer-based Lanchester-type models of 
quite complex combat systems have been developed. Such models currently exist for 
almost the entire spectrum of combat operations, from combat between battalion-sized 
[9] and division-sized [10] units to theater-level operations [12,14]. 

A simple combat model, however, may yield a clearer understanding of important 
interrelationships that are difficult to perceive in a more complex model, and such 
insights can provide valuable guidance for more detailed computerized investigations 

(see [7,28]). In this paper we present a new important numerical procedure that 

§ 

facilitates parametric analysis of battle outcomes for such simplified Lanchester- 
type models of combat between two homogeneous forces with temporal variations in each 
side’s fire effectiveness. Previously, such battle-outcome information could only be 
readily obtained from constant-coefficient models These results are not only 
significant in their own right but are also useful in the quantitative analysis of 
time-sequential combat strategies (see, for example, [19-20]). 

It is important for the military operations analyst to have a clear understand- 
ing of how force-level and weapon-system-performance parameters interact to determine 

^Corresponding stochastic formulations are for all practical purposes analytically 
intractable (see Note 1 on p. 65 of [22]). 

^Also frequently called differential models of combat [10] . 

§ 

In particular, the parametric examination of force-annihilation prediction. 

Bonder [1-3] has emphasized the deficiencies of constant-coefficient models 
(see Section 3 below) . 
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a battle’s outcome. Such knowledge is particularly useful in weapon-system and force- 
level planning activities for defense planning.^ S. Bonder’s [1,3,7] pioneering work 
on methodology for the evaluation of military systems (particularly mobile systems 
such as tanks, mechanized infantry combat vehicles, etc.) provides a motivation for 
interest in variable-coefficient, deterministic, Lanches ter- type combat models such 
as we consider in this paper. He has stressed (see pp. 30-31 of [7]) the importance 
of analytical solutions to such models for developing insights into the dynamics of 
combat by portraying the relation between various factors in the combat attrition 
process and the surviving numbers of forces and for facilitating sensitivity and 
other parametric analyses (see [5]). Unfortunately, as work by Bonder and Farrell 
[7] and Taylor [18,22] shows, the analytical (i.e. infinite series) solution to 
variable-coefficient equations generally by itself^ provides little information 
about battle outcome because of its complexity. Therefore, one must seek new ways 
for developing insights . ■ 

Taylor and Comstock [24] have given results that allow one to predict battle 
outcome (in particular, force annihilation ) in theory without having to spend the time 
and effort of computing force-level trajectories. To be computationally practical, 
however, their results require the determination of the so-called parity-condition 
parameter (”the enemy force equivalent of a friendly force of unit strength'*), which 
depends on only the model’s attrition-rate coefficients. They analytically determine 
the parity-condition parameter for power attrition-rate coefficients with "no offset," 
which allow one to model combat between two weapon systems with the same maximum 

^Especially since one frequently uses models that are so complicated that trends are 
not directly discernible without extensive (and costly) computer runs. 

^I.e. without explicitly computing force-level trajectories. 

^Bonder and Honig [8] point out, however, that force annihilation may not always 
be the best criterion for evaluating military operations. See pp. 192-242 of Bonder 
and Farrell [7] for a detailed Lanchester-type analysis of an attack scenario for 
which other ’’end of battle conditions" play the principal role. Nevertheless, it 
is of considerable interest (especially for developing insights into the dynamics 
of combat) to be able to easily predict the occurrence of force annihilation. 
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effective range but different range dependencies for each system* s fire effective- 
ness (see also [23]). It is the purpose of this paper to show how to determine the 
parity-condition parameter in other cases, in particular for power attrition-rate 
coefficients with ’’positive offset,” which allow one to model such combat between 
weapon systems with different maximum effective ranges. Our results allow one to 
study in general such variable-coefficient combat models almost as easily and 
thoroughly as Lanchester’s classic constant-coefficient model. 

The organization of this paper is as follows. We first review Lanchester-type 
equations of modern warfare, especially variable-coefficient formulations. Next we 
review force-annihilation-prediction conditions for such models and show how to use 
knowledge about the parity-condition parameter for one set of attrition-rate 
coefficients to numerically determine it in related cases of interest. This general 
theory is then applied to the important case of offset power attrition-rate 
coefficients, with detailed numerical examples being given. 

2. Lanchester’s Classic Combat Formulation . 

F. W. Lanchester [16] hypothesized in 1914 that combat between two military 
forces could be modelled by^ 

dx/dt = - ay , dy/dt = - bx , (1) 

with initial conditions 

x(t = 0) = x^ , y(t = 0) = y^ , (2) 

where t = 0 denotes the time at which the battle begins, x(t) and y(t) denote 
the numbers of X and Y at time t, and a and b are nonnegative constants which 
are today called Lanchester attrition- rate coefficients and represent each side’s 
fire effectiveness. We will refer to (1) as Lanchester ’s equations of modern 

^See also p. 45 of [22]. 

^The equations are only valid for x, y > 0. For example, the first becomes 
dx/dt = 0 when x = 0. 
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warfare.^ Various sets of physical circumstances have been hypothesized to yield 
them: for example, (a) both sides are aimed fire and target acquisition times are 

constant [27], or (b) both sides use area fire and a constant density defense (see 



p. 345 of [13]). 



From (1) Lanchester deduced his famous square law 



b(xQ - x^(t)) = a(yQ - y^(t)) . 



(3) 



Consider now a battle terminated"^ by either force level reaching a given "breakpoint" 



for example, Y wins when - x(t^) - x^p 



BP 

^0 ^f ^ ^BP 



rBP u 

fy ^0’ 



tf, x^, y^ denote final values and x^p denotes X’s breakpoint which is a given 

fraction f®^ of his initial strength. It follows from (3) that 
X 



Y will win if and only if 

^0 



/aU - (fy®)^} 
bU - (ff )2} 



(4) 



BP BP 

which for a fight-to-the-finish (i.e. f^^ “ ^y “ becomes the classic result 



Y will win a fight-to-the-finish if and only if 






(5) 



Unfortunately, no relationship similar to (3) holds in general for variable 
attrition-rate coefficients.® This paper, nevertheless, shows how (5) generalizes 
in these cases. Recalling that the time history of the X force level is given by 

x(t) = Xq coshv^ t - y^ sinhv^^ t , (6) 

we see that the battle trajectories depend on the two weapon- system-performance 



^Other forms of Lanchester- type equations appear in the literature [13,18], but we 
will not consider these here. 

^The modelling of battle termination is a problem area in contemporary defense 
planning studies (see pp. 524-525 of [25]). 

§ 

Except when a(t)/b(t) = constant (see p. 48 of [22]). 

II So far we have not been able to generalize (4). 
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parameters: (I) the intensity of combat v^, and (II) the relative fire effectiveness 

a/b. Only the relative fire effectiveness, however, determines the battle’s outcome 
[see (4) and (5) above]. 



3. Variable Attrition-Rate Coefficients 

Bonder [2] has pointed out that in many cases (for example, in the case of mobile 
weapon systems) the validity of the assumption of constant attrition-rate coefficients 
is open to question (see also [1,3,7]). Thus, we consider 

dx/dt = - a(t)y, dy/dt = - b(t)x , (7) 

where a(t) and b(t) denote time-dependent attrition-rate coefficients. We assume 
that a(t) and b(t) are defined, positive, and continuous for t < + <» with 

^ also assume^ that a(t), b(t) e L(tQ,T) for any finite T. We further 

take a(t) and b(t) to be given in the form a(t) = k^g(t) , b(t) = k^h(t) , where 
k^, k^ are positive constants chosen so that a(t)/b(t) = when g(t) = h(t) . 

Analogous to the constant-coefficient case [see discussion after (6)], we have the 
two weapon-system-performance parameters: (I) the intensity of combat, I(t) = 

/a(t) b(t) ; and (II) the relative fire effectiveness, R(t) = a(t)/b(t). We accord- 
ingly introduce the combat-intensity parameter and the relative-fire-effective- 

ness parameter X^ defined by 



and 



X = k /Ic 
R a b 



( 8 ) 



Two significant developments in the Lanchester theory of combat during the 
1960’s were the development of methodology for (a) the prediction of Lanchester 
attrition-rate coefficients from weapon-system-performance data by S. Bonder [2,4], 
and (b) the (maximum likelihood) estimation of such coefficients from Monte Carlo 
simulation output by G. Clark [11]. Both these developments and others^ have 



t. 



For convenience, we introduce the notation that a(t) C L(t^,T) means that 
T 

/ a(t)dt exists (and is given by a finite quantity) . From our assumptions about 
a(t) and b(t), it follows that a(t) i L(tQ,T) implies that a(t)dt = + “ . 

See [22] for further references. 
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generated interest in the model (7) and facilitated its application (and that of 
its generalization to combat between heterogeneous forces [7]) to defense planning 
studies . 

A large class of tactical situations of interest can be modelled with the 
following general power attrition-rate coefficients [7,22,24] 

a(t) = k^(t + C)^, and b(t) = k^(t + C + A)'^, (9) 

where A, C ^ 0. We will call A the offset parameter , since it allows us to model 
(with y, V ^ 0) battles between weapon systems with different maximum effective 
ranges. We will call C the starting parameter , since it allows us to model (again 
with y, V battles that begin within the minimum of the maximum effective ranges 

of the two systems. For example, let us consider Bonder’s [1,3] model of a constant- 
speed attack on a static defensive position (see also [18,22]). Then we have 

dx/dt = - a(r)y , dy/dt = - 3(r)x , (10) 



where r(t) = - vt denotes the distance (range) between the two opposing forces, 

Rq denotes the battle’s opening range, v > 0 denotes the constant attack speed. 



ot(r) = 



o.o(l-r/R^)>' 



for r > R , 
— a' 



for 0 < r < R , 
— — a’ 



( 11 ) 



y ^ 0, and R^ denotes the maximum effective range of Y’s weapon system. Similarly 
for B(r), with exponent V ^ 0. In (11) the parameter y allows us to model the 
range dependence of Y’s fire effectiveness (see Fig. 1). The offset and starting 
parameters are given by 



A = (Rg - R^)/v , and C = (R^ - Rq)/v , (12) 

and the assumption A, C > 0 implies that R„ > R > Rr,» From considering (12) 

— p — ot — u 

and Fig. 2, the reader should have no trouble understanding our terminology for A 
and C . 

The time history of the X force level, i.e. the solution x(t) to (7), is 



given by [22] 
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Rotf 




of battle is denoted as = 1250 meters and (as shown) 
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(13) 



x(t) = XQ{CY(0)Cj^(t)-SY(0)S^(t)} - {Cj^(0)S^(t)-S^(0)C^(t) } , 

where the hyperbolic-like general Lanchester functions (GLF) C^,(t) and S„(t) are 
linearly independent solutions to the X force-level equation 




with initial conditions 



1 da ) dx 
a(t) dt / dt 



a(t) b(t)x = 0 , 



(lA) 



C^(t tg) 1 , 



Sx(t t^) 0 , 



(15) 



{[l/a(t)]dC /dt(t)} = 0 , {[l/a(t)]dS /dt(t)} = 1/v^ , 

X t-tp X t-tQ R 

where t^ denotes the largest finite time at which a(t) or b(t) ceases to be 
defined, positive, or continuous. For example, t^ = -C for the general power attti- 
tion-rate coefficients (9). The time history of the Y force level may be similarly 
obtained, with C^Ct) and S^Ct) being analogously defined for the corresponding 
Y force-level equation. 

For the numerical determination of the parity-condition parameter, it is con- 
venient to introduce a new independent variable s defined by 



t 

s = KA^ g(o) do , 



(16) 



where the parameter K is to be chosen to simplify the form of J(s) given by (18). 
We denote s(t = 0) as s^, and then ^ 0 if and only if 0. The sub- 

stitution (16) transforms (lA) into the normal form [17] 



where 



^ - J(s)x = 0 , (17) 

ds 




( 18 ) 



and t = t(s). We also define the normal-form hyperbolic-like GLF c^(s) and 

s (s), which satisfy (17) and the initial conditions 
A 
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c^(s = 0) = 1 



Sx(s = 0) = 0 , 



( 19 ) 



dc^/dsCs = 0) = 0 , 

It follows that 

c^(s(t)) = C^Ct) , and 



ds^/dsCs = 0) = 1 . 
Sx(s(t)) = KS^Ct) . 



( 20 ) 



4. Force-Annihilation-Prediction Conditions 

THEOREM 1 (Taylor and Comstock [24]): Assume that either a(t) it L(0,+ 

b(t) it L(0,+ “) . Then the X force will be annihilated in finite time 



Xq ^ ( C (0) - Q*S (0) 

r" H —\- 

^0 ( Q Cy(0)- 8^(0) 



oo) or^ 

if and only if 
( 21 ) 



where the parity condition parameter Q unique and given by 



■ g* 



Remark 1 : We also have lim {S (t)/Cy(t)} = Q 



( 22 ) 



Remark 2 : When t^ = 0, (21) simplifies: X will be annihilated in finite time if 

and only if x^/y^ < »^/Q • 

Remark 3 : The result (22) suggests a numerical procedure for approximately deter- 

mining the parity-condition parameter Q : we may approximate the parity-condition 

parameter Q* by Q = l/{S^(t) /C^(t) } , where t is a "suitably large" value of t. 
In other words, we may estimate Q* simply by picking a large value for t (we 
denote this selected large value by t), computing S^(t) and C^(t), and then form- 
ing their ratio. Our estimate for Q is then given by Q = 1/{S^( t) /C^( t) } . The 

/\ 

only problem is that we don’t know how large to take t for "satisfactory" estima- 

‘ "k /V 5/lf /S /N 

tion of Q : There is an estimation error, E(t) = Q - Q(t), which depends mono- 
tonically on t, and a priori we don’t know how large this error is. The present 



paper develops a bound on the magnitude of this error, and our new error estimate 
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allows the goodness of approximation to be easily evaluated in many cases of interest. 

We may also determine the parity-condition parameter with the normal-form hyper- 
bolic-like GLF, since lim {s (s)/c (s)} = 1/Z = K/Q , where Z is called the 

S”^ ~r A A 

modified parity-condition parameter . In fact, we will find it more convenient to do 
so. With this in mind, let us introduce the Y-functions c^(s) and s^(s) [corres- 
ponding to c^(s) and s^(s)] defined by 

dc^/ds = J(s)s^ , ds^/ds = J(s)c^ , (23) 

with initial conditions 

c^Cs = 0) = 1 , s^Cs = 0) = 0 . (24) 



It follows that 
Y equation 

and 



and 3 ^( 5 ) are linearly independent solutions to the modified 



A. i 

ds (J(s) ds / 



- y = 0 , 



(25) 



CY(s(t)) = C^Ct), s^CsCt)) = (1/K) S^Ct) . 



(26) 



In terms of the new time variable s defined by (16), Theorem 1 reads as follows 



THEOREM 2: Assume that either a(t) t L(0, + «>) or b(t) t L(0, + <») . Then the X 



force will be annihilated in finite time if and only if 



^0 ^ 1 Z Cy(Sq) - Sy(Sq) 



(27) 



where the modified time variable s is given by (16), and c (s) , s (s) , c (s) , and 

A A 1 

Sy(s) denote the normal-f orm hyperbolic-like GLF . The modified parity-condition 
parameter Z unique and given by 



lim 

S ^ + 00 



Cj^(s) 




(28) 



We observe that . . 

•k k 

Q = KZ , 



(29) 
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and lim , {s„(s) /cv(s) } = Z . When (27) holds, the time to annihilate X, denoted 

S "TOO I 1 ' 

X X 

as t , is determined by x(t ) = 0. If we denote the quotient of the two normal- 
a a 

form hyperbolic-like GLF s^(x) and c (s) as n„(s) , then it follows from (13) that 

A y X 



rix(s(t )) = — , 

{xqSy(So) + yo(A:^/K)Cjj(sQ)} 

where 



njj(s) = Sj^(s)/Cj^(s) . 



(30) 



(31) 



5. Determination of the Parity-Condition Parameter . 

* 

We will now show how knowledge about the modified parity-condition parameter Z 
for one pair of attrition-rate coefficients, a(t) and b^(t), allows us to determine 
Z for a related pair, a(t) and b(t). With this in mind, let us denote c (s) 

A 

corresponding to a(t) and b(t) as c (s;a,b), and similarly for s and r] . 

X XX 

In other words, we will now write (31) corresponding to the attrition-rate coefficients 
a(t) and b(t) as 

nj^(s;a,b) = Sj^(s;a,b)/Cj^(s;a,b) . (32) 



In this notation, we will write (28) as 



t 



lim 
s ->■ +" 



Tix(s;a,b) 



1/Z*[a,b] . 



(33) 



Our main result is Theorem 5, which gives an error estimate for the approximation 

* 

that we propose for Z . The theoretical basis for Theorem 5 is given by Theorem 4, 
which (in turn) is a consequence of Theorem 3. The proof of Theorem 3 follows along 
the lines of well-known arguments (see p. 225 of [15]). 



"t" ^ r 

We use the notation Z [a,b] to show that the modified parity-condition parameter is 
a functional (i.e. a function for which the independent variables themselves are 
functions), which depends on only the attrition-rate coefficients a(t) and b(t). 

In other words, the attrition-rate coefficients are functions defined for tQ £ t < +“ , 
and the parity-condition parameter depends on these entire functions (and not 
merely particular values of them) . 
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THEOREM 3 (Comparison Theorem): Let x(t) and x^(t) satisfy 






1 



dt i a(t) dt 



bj^(t)x^ = 0 , 



with initial conditions 



X ( t tQ ) CL f 



Xi(t = tp) = a , 



Ul/a(t)]dx/dt(t)} = 3 , 

0 



{ [l/a(t) ]dx /dt(t)} = 3 , 

0 



where a(t) > 0 and b^(t) < b(t) for all t > t^. Then x^(t) < x(t) for all 
^ ^ as long as x(t) > 0. 

* 

The basic theoretical result upon which our numerical determination of Z 
is based is 

THEOREM 4. Assume that b^(t) < b(t) for all t > t^. Then 

nj^(s;a,b) < 1/Z*[a,b] < rix(s;a,b) + { (1/Z*[a,b^ ] ) -rij^(s ;a ,b^) } . (3A) 

PROOF: We observe that [24] n^^(s;a,b) satisfies the differential equation 

A 



drij^/ds(s;a,b) = l/{c^(s ;a ,b) } , 



(35) 



with r)v(s ” 0;a,b) = 0, and similarly for r)v(s ; a ,b-, ) . Theorem 3 (the comparison 
X XI 

theorem) yields that c (s;a,b) > c^(s;a,b ) for all s > 0. Thus, for all s > 0 

X XI 

drij^/ds(s;a,b) < drij^/ds(s;a,b^) , 



whence integration between 0 and s yields the desired result. 



Q.E.D. 



Similar to the observations made in Remark 3 above, we observe that (33) 
suggests that we estimate Z [a,b] with Z defined by 

Z(s;a,b) = l/rij^(s;a,b) , (36) 

/\ 

where s denotes a suitably chosen value for s. Moreover, from (35) we see that 
r|^(s;a,b) is a strictly increasing function of s so that the larger we take s 
in (36), the better our approximation becomes. The only problem (see Remark 3) is 
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that a priori we don’t know how large to take s for ’’satisfactory" estimation 

* /s 

of Z . Theorem 5, however, tells us exactly how large to take s. 

THEOREM 5 (Error Estimate for Approximation): Assume that b^(t) < b(t) for all 

* 

t > t^. Let denote the fractional error made in the estimation of Z [a,b] 

by Z(s;a,b), i>e > 

b) - z* l a,bl ^3,^ 

® 2 [a,b] 

Then 

0 < fg(s) < {(1/Z*[a,b^]) - Dx(s;a,b^)} * Z(^,a,b) . (38) 

ys 

PROOF: The theorem follows by simple algebraic manipulation after setting s = s 

in (34) and using (37). Q . E.D. 



Thus, we have presented a method for numerically determining Z [a,b]. We 
simply pick a large value for s (we denote the selected value as s) , compute 

/V /S /N ^ 

s (s) and c (s) , and then compute the estimate Z(s;a,b) according to (36). 

\ A. 

Theorem 5 allows us to know the accuracy of our approximation, which can be improved 

^ * 
by taking s larger. Thus, we can numerically determine Z [a,b] to any specified 

* 

degree of accuracy once Z [a,b^] is known. In the next section we apply this 
theory to the analysis of battles modelled with offset power attrition-rate 
coefficients . 



6. Application of Theory to Offset Power Attrition-Rate Coefficients 

In the application of Theorems 4 and 5, two pairs of attrition-rate coefficients 
are involved: one pair for which the modified parity-condition parameter is known 

[denoted as a(t) and b^(t)], and one for which it is to be determined [denoted 
as a(t) and b(t)]. Accordingly, we rewrite (9) with A > 0 as 

a(t) = k^(t + C)^, and b(t) = k^(t + C + A)^, (39) 



12 



where (as before) C ^ 0. We will refer to these coefficients (39) for which A > 0 
attrition- rate coefficients with ’’ positive offset .” If we choose 

K = [Xj/(u+l)]^P"^ , (40) 

it follows from (16) that the modified time variable s is given by 

S(t) = [X^/(y+l)]^P(t + 0^“^^ , (41) 

and the invariant J(s) of the normal form (17) simplifies to 

J(s;a,b) = J(s;Y,y,v) = s^(l + ^) , (42) 

s 

? IT 

where p = (v+l)/Z, a = l/(y+l), B = (v-y)/(y+l), y = A-[X^/(y+l)] , and Z = y+v+ 2. 

Here we have denoted the invariant corresponding to the attrition-rate coefficients 
a(t) and b(t) as J(s;y,p,v), since we may take y, y, and v as a basis for 
generating the four parameters 3, y, and V that explicitly appear in the right- 

hand side of (42). Furthermore, we will denote the normal-form hyperbolic-like GLF 
that correspond to J(s;y,y,v) as c^(s;y,y,v) and s^(s;y,y,v). 

The known results [24] that we use in the Theorems 4 and 5 are for the case 
of power attrition- rate coefficients with no offset [i.e. set A = 0 in (9)] 

a(t) = k^(t + C)^, and ^l^^^ ” 

where C ^ 0. We observe that b^(t) < b(t) for all t > -C. It follows that 
J(s;a,b^) = and [24] 

Z*[a,bJ = r(l-p)/r(p) . (44) 

"k k 

Thus, for the bound on Z [a,b] = Z (y,y,v) given in Theorem 4 and the error estimate 
for our approximation (36) given in Theorem 5, we have [23] 

rij^(s;a,bj^) = p^^"^P^ T^(S) , (45) 

where S = 2p q = 1-p, and T denotes a Lanchester-Clif ford-Schafli 

i 

(LCS) function, which is analogous to the hyperbolic tangent (see Table I) . 

^These functions were introduced in [22] and redefined for reasons of force-annihi- 
lation prediction in [23]. 
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TABLE I. 



LANCHESTER-CLIFFORD-SCHAFLI FUNCTIONS 



“ (x/2)2>^ 



Ft,(x) - r(o) 1(1 r(i( + 



“a'*) - '■(“> W T 1) 



2(k+a) 



T (x) = H, (x)/F (x) 
a' 1-a' 'a 



Relation to normal-form GLF . 

c (s) = F (S) 

X q 

Cy<s) = Fp(S) 



s^is) = 
8^(8) = 



where q = l~p and 

S(s) = 2ps^^<^P> 



NOTE: For y = V, we have 



(I) 


C^(s) = 0^(8) 




(II) 


S^(S) = 8^(8) 




(III) 


\(s) = n^(s) 


= Tj^/ 2(8) = tanh s. 



In fact, ‘ 



n(s;a,b^) = tanh s, when y = v . (46) 

We have thus shown that the following theorem holds. 

THEOREM 6: Assume that either y > -1 V > -1. Then for a battle modelled with 

the offset power attrition-rate coefficients ( 39 ) , bounds on the modified parity - 
condition parameter Z (y,y,v) are given for y ^ 0 by 



n^(s;Y,y,v) < 



z (Y.y.v) 

l/(2p) 



< rix(s;Y.y,v) + p 



q-p 



I -’■<,«>} ■ 



(A7) 



where q = 1-p, S = 2ps , and r)y(s;y,y,V) denotes the quotient of two normal - 



form hyperbolic-like GLF for the attrition-rate coefficients (39), i.e. ri„(s;y,y,v) 

A 

= s^(s;Y,y,v)/c^(s;Y,y,v) . 



It follows from Theorem 6 (or, equivalently, Theorem 4) that if we approximate 
Z*(Y,y,V) with Z(s;Y,y,V) defined by 

Z(s;Y,y,v) = l/Dj^(s;Y,y,v) , (48) 

then bounds on the fractional error made in this estimate are given by 

0 < fg(s) < T^(S)| nx(s;Y,y,v), (49) 

where f^(s) denotes the fractional error and is defined by (37) • 

The right-hand inequality in (49) [equivalently, (47)] tells us exactly how 
large to take s for the estimation of Z (y > 0,y,v) by Z(s;y,y,v) to any 
specified degree of accuracy. The LCS function T^ is involved in the bound on 
the fractional error 1,^(3) in this estimate when y v.^ Thus, the LCS functions 

£i 

as redefined by Taylor and Brown [23] yield valuable information about battles 
modelled with not only the power attrition-rate coefficients with no offset (43) but 
also the offset power attrition-rate coefficients (39) . Availability of tabulations 
I of these LCS functions is discussed in [23]. 

^This result is one of our reasons for introducing the normal form (7). 

^As (46) and Table I show, T^(S) = tanh s when y - V. 
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7. Numerical Results 



In this section we will examine a couple of numerical examples to show how 

* 

the modified parity-condition parameter Z may be numerically determined and to show 
some important insights into the dynamics of combat that may be consequently obtained. 
In order to numerically determine the modified parity-condition parameter for the off- 
set power attrition-rate coefficients (39), we must use knowledge about how quickly 
the limiting value (i.e. Z [a,bj^]) of a hyperbolic- tangent-like function of a related 
pair of power attrition-rate coefficients with "no offset" (43) is reached as its 
argument is increased [recall Theorem 6 and (49)]. In Fig. 3 we see that this limit- 
ing value, denoted as Z (y,v) = Z [a,b^], is quite quickly reached: if one takes 

s = 10.0, then Z (y,v) is approximated to better than six decimal places by 

/N A . ^ 

Z(s;y,v) = l/r|^(s;y , v) , where is given by (45). Experimental computing for 

* 

various values of y and V and comparison with the known value (44) for Z (y,v) 
bears out this degree of accuracy [i.e. speed of convergence of Z(s;y,v) to Z ] 

A A 

for essentially all allowable values of y and V. Thus, Z(s;y,v) for the coef- 

* A 

ficients (43) has essentially converged to Z (y,v) when s = 10.0, and by Theorem 6 

/S A 

or (49) we know that the same is true for Z(s;YjyjV) for the coefficients (39). 

We have accordingly generated by this procedure the results shown in Fig. 4. 
For computing f)^ = s^/c^, we have used the series solutions shown in Tables II and 
III. [In Tables II and III we have for convenience denoted, for example, s (s;Y>y>v) 

A. 

simply as s^(s;y,v), i.e. s^(s;y,v) denotes s^ corresponding to the general 
power attrition-rate coefficients (9) with exponents y and v.] The series were 
obtained by solving (17) by the method of successive approximations (see [18]). We 
used these series instead of developing approximate solutions by finite-difference 
methods because we did not have any error bounds for the latter. 

Let us now give an intuitive interpretation of the curves shown in Fig. 4 

* 

of the modified parity-condition parameter Z plotted versus the modified offset 

* 

parameter Y. In Taylor and Comstock [24] it is shown that Z may be considered 
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TABLE II 



NORMAL- FORM OFFSET LINEAR LANCHESTER FUNCTIONS 



Cj^(s;l,l) = I 



00 2k k . / \ 1 

s ^ .J JLJ 



k£o (2k): \ ^ ^3 



X 



00 2k+l k 
(s;l,l) - ^ / 9 V 4 . 1 1 ' I 






where the offset coefficients are given by 






1, 



and for k > 1 



4 



4k(4k-2) 

(4k-j)(4k-2-j) 





for 0 _< j _< k 



( 

I 



0 



1, and for k ^ 1 



4k(4k+2) 

(4k-j)(4k+2-j) 



{b 



j 

k-1 






for 0 £ j £ k 



We have adopted the convention that 



A^ 

k 




for j < 



0 or j > k. 



TABLE III 



OFFSET POWER LANCHESTER FUNCTIONS FOR y = 1 AND v = 2 



'X' 



c„(s;l,2) = r(|) I ^ 



,2 5/4, 2k 

V c s ) 



k=0 kl T(k 



.h j-o 



X 



c (2 5/4 2(k+2/5) 

2^1/5 r.2, V ^ f bJ 



(s;l,2) = (1)^''^ r(|) I 



1^=0 k.' r(k + j) j=o 



Y 

7T 



where the offset coefficients are given by 



t 



= 1, and for k >. 1 



J = 5k(5k-2) . j-2 j-1 j , for 0 < -i < 

\ (5k-j)(5k-2-j) ^^-1 ^-1 ^-1^ u < J < 



2k 



= 1, and for k > 1 



B? = 



k (5k 



-IhS'Lj) t'klj + bJ.^} for 0<3<2k 



^We have adopted the convention that A^, B^ = 0 for j < 0 or j 



> 2k. 



(50) 



to be the initial Y force level that leads to a draw^(i.e. parity between the 

forces) in the battle against an X force of “unit strength" 

dx/ds = -y with x(s = 0) = 1 , 

dy/ds = -J(s)x with y(s = 0) = Z* , 

where J(s) denotes the invariant of the normal form (17). Thus, we may consider 

Z to be “the Y equivalent of an X force of unit strength" for the modified 

battle (50). Now let us consider the general power attrition-rate coefficients (9) 

with exponents y and V. As we did in Tables II and III, we will denote the corre- 

* * 

spending J as J(s;y,v) and Z as Z (y,v) to stress the dependence on y and 

V (but suppressing that on y). We then have from (42) that J(s;l,l) = 1 + y//s 

and J(s;l,2) = /s (1 + y/v^)^. From (44) we find that Z (1,1) = 1.000 and Z (1,2) 

= 0.806 for y = 0. Observing that for y > 1 we have J(s;l,l) < J(s;l,2) for all 

* 

s ^ 0, it is intuitively clear from (50) and the interpretation of Z as a force 

"k k 

equivalent that we must have Z (1,1) < Z (1,2) for all y > 1 because X always 

has greater fire effectiveness against Y when y = 1 and V = 2 than when y = 1 

and V = 1. However, for y near zero, the situation is reversed and Z (1,2) must 
lie below Z (1,1) for y near zero. Thus, we have given an intuitive explanation 

* k 

of why Z (1,2) lies below Z (1,1) for y near zero but above it for y > 1 as 

Fig. 4 shows. 

Next, we will consider numerical results for a particular battle to show some 
of the important insights that may be gained into the dynamics of combat from our new 
results. As in [18,22,23] we consider S. Bonder’s [1,3] model (10) for the constant- 
speed attack of mobile forces against a static defensive position. We will focus on 
the new results given in this paper (in particular, the prediction of battle outcome 
from initial conditions without explicitly computing the force-level trajectories). 
Input data and computed parameter values are shown in Table IV. We will now consider 
two cases: (I) = 1500 meters, and (II) Rq = 1250 meters. 



In other words, x(s) and y(s) > 0 for all s e [0,+oo) but lim^_^_^x(s) = 0 
= lim . y(s). 
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TABLE IV 



PARTICULARS FOR THE NUMERICAL EXAMPLES 



1. Input Data 
y = V = 1 

a^= 0.06 X casualties/minute/Y unit 

3q= 0.6 Y casualties/minute/X unit 

R = 1500 meters, Rq= 2000 meters 

a 3 

V = 5 miles/hour 

2. Parameter Values 

-3 

k = 5.364 X 10 X casualties/minute/Y unit 
a 

-3 

k^= 4.023 X 10 Y casualties/minute/X unit 
P = q = 1/2 

1/2 

A = 3.728 minutes, Y = 0.320 (casualties ’minutes) 



When Rq = 1500 meters, we have C = 0 and = 0. The maximum time that 

the battle can last is t = 11.18 minutes, since at this time the advancing 

max ^ 

"k ic 

attackers overrun the defensive position. In this case Z (y,y,v) = Z (0.32,1,1) 

= 1.381, so that Theorem 2 tells us that X can be annihilated <==> ^q/Yq < 0.264. 

X 

By (30) the X-force annihilation time is given by )) “ 2.739x^/y^. For 

A a u u 

X f 

Xq = 10 and y^ = 50, we have n^(s^) = 0.54772 so that by the techniques introduced 

in [23] we find s^ = 0.771. Hence, (36) yields t^ = 10.25 minutes and r^ = 125.7 
a a a 

meters. Further results are given in Table V. 

When Rq = 1250 (see Fig. 5 of [22]), we have C = 1.864 minutes, s^ = 0.0255 

and t =9.32 minutes. In this case X can be annihilated <==> x^/y^ < 0.281 
max 0-^0 

with the X-force annihilation time given by q (s^) = (l.OOlu^. + 0 .009) / (0 . 127u + 0 . 366) , 

A a U u 

where uq = ^o^^O* Numerical results are given in Table VI. Finally, these parametric 
results should be contrasted to those previously possible (e.g. compare them with, 
for example, the single force-level trajectory for Rg = 2000 meters shown in Fig. 5 
of [22]). 

8. Discussion 

S. Bonder [1-3] has emphasized the shortcomings of constant-coefficient 
Lanchester-type combat models. Work by Bonder [1,2], Clark [11], and others [7] on 
the prediction of Lanchester attrition-rate coefficients^ has generated interest in 
variable-coefficient models. Moreover, there is not only intrinsic interest (see 
[1,3]) in the model (7) but also interest for obtaining insights into the behavior 
of complex Lanchester-type system models§ that have been enriched in military detail 

^These computations involve the generation of a table of s , c^, and Hy for 
y = 0.32, y = V = 1 (see [23]). 

^See Taylor and Brown [22] for further discussion and references. 

§For example, the Bonder-IUA model (see [7-9]). 
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TABLE V 



ANNIHILATION OF THE X FORCE AS A FUNCTION 
OF THE INITIAL FORCE RATIO FOR R^ = 1500 METERS 



(xo/Yo) 


X 

(minutes) 


(meters) 


0.250 


14.09 


t 


0.200 


10.25 


125.7 


0.167 


8.80 


319.4 


^t =11.18 minutes and 

max 


x^ = x(r = 0) = 2.48 





TABLE VI 

ANNIHILATION OF THE X FORCE AS A FUNCTION 



' THE INITIAL 


FORCE RATIO FOR 


R =1250 METERS 
0 


o 

o 


X 

t^ (minutes) 


r^(meters) 


0.250 


10.87 


t 


0.200 


8.17 


154.4 


0.167 


6.93 


320.4 



=9.32 minutes and = x(r = 0) =1.74 
max f ' ' 



see [7-10,12,14]). The attrition-rate coefficients in (7) represent the fire 
effectiveness of the combatants and allow us to model temporal variations in fire 
effectiveness on the battlefield. Interest in the general power attrition-rate 
coefficients (9) is provided by S. Bonder’s [1,3,8] constant-speed attack model^ 
(lO)-(ll) and his examination of the range dependence of attrition-rate coefficients 
for various weapon systems (see pp. 196-200 of [7]). 

We have given results that allow one to study the variable-coefficient model 
(7) [especially with the general power attrition-rate coefficients (9)] almost as 
easily and thoroughly as Lanchester’s classic constant-coefficient model (1). Taylor 
and Comstock [24] (see Theorems 1 and 2 above) have shown how to predict force 
annihilation without having to spend the time and effort of explicitly computing 
force-level trajectories. Using their theoretical results, we gave results in a 
previous paper [23] that made combat modelled by power attrition-rate coefficients 
with no offset [i.e. A = 0 in (9)] almost as easy to analyze as the constant- 
coefficient case. The results of the paper at hand allow one to analyzed combat 
modelled by power attrition-rate coefficients with positive offset [i.e. A > 0 
in (9)] just as conveniently. 

Theorem 1 (see also Theorem 2) is the generalization of the classic constant- 
coefficient result (5) to cases of time-dependent attrition-rate coefficients. How- 

■k 

ever, one needs to know the value of the so-called parity-condition parameter Q 
in order to predict force annihilation in specific instances. In this paper we have 
presented theoretical considerations (see Section 5 above) for the noniterative 
numerical determination of the parity-condition parameter. We applied our general 
theory to the specific case of general power attrition-rate coefficients (9) (see 

^Thus, the range between firer and target changes during the engagement. 

^Modellipg, for example, combat between two weapon systems with the same maximum 
effective range. 

§ 

Modelling, for example, combat between two weapon systems with different maximum 
effective ranges. 



18 



Section 6) and illustrated these theoretical results by considering some numerical 
examples (see Section 7) . 

Curves of the modified parity-condition parameter Z plotted against the 
modified offset parameter y such as those shown in Fig. A allow one to parametrically 
analyze "modern" combat modelled with the general power attrition-rate coefficients 
(9). For example, we can now parametrically^ determine whether the defender will 
be overrun in Bonder's [1,3,8] constant-speed-attack model (10) with attrition-rate 
coefficients (11) without having to compute the entire force-level trajectories. We 
illustrated this analysis capability with some numerical examples, which showed that 
the defender's annihilation (i.e. saturation of his defensive position with offensive 
fire) depended on the initial force ratio (of defender to attacker) being below a 
certain threshold value. Our new results allow one to not only easily determine such 
force-ratio thresholds of survivability but also study their dependence on weapon- 
system-capability parameters. 

Our new results let us conveniently obtain much valuable information about 

the model (7).^ Previously one was limited to only computing force-level trajectories, 

but now we can predict battle outcome (in particular, force annihilation) without 

explicitly computing force-level trajectories. Moreover, these new results facilitate 

§ 

parametric analysis of such combat situations. In particular. Theorems 1 and 2 
explicitly exhibit a tradeoff between quality (as quantified by the relative-f ire- 

ic 

effectiveness parameter and the parity-condition parameter Q ) and quantity 

R 

(as quantified by the initial force ratio ^q/Yq) weapon systems in combat 

against each other. In other words, one can use an expression like (21) to develop 

^Varying, for example, the maximum effective range of the defender's weapons. 

^The classic ordinary differential equation theories (see, for example, Ince [15]) 
were inadequate to answer many important questions (for example, "Who will win? 

Be annihilated?") about such combat models. 

§ 

S. Bonder [5] has suggested that an increased emphasis be placed on parametric 
analyses in systems analysis studies (see pp. 21-22 of [5]). 
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quantitative insights into how the quality of a weapon system may be substituted 

for sheer numbers. Moreover, an unanswered theoretical question is to determine 

* 

how the parity-condition parameter Q depends on the combat-intensity parameter 

Xt Sind the relative-fire-effectiveness parameter . Finally, our results here 
I ^ 

are signposts as to the difficulty of analytically extracting information (particu- 
larly parametric information without excessive computations) from variable-coef- 
ficient Lanchester-type models such as (7). 
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